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Abstract: We derive a set of first-order differential equations obeyed by the S-matrix of 
planar maximally supersymmetric Yang-Mills theory. The equations, based on the Yangian 
symmetry of the theory, involve only finite and regulator-independent quantities and uniquely 
determine the all-loop S-matrix. When expanded in powers of the coupling they give deriva- 
tives of amplitudes as single integrals over lower-loop, higher-point amplitudes/ Wilson loops. 
We outline a derivation for the equations using the Operator Product Expansion for Wilson 
loops. We apply them on a few examples at two- and three-loops, reproducing a recent result 
on the two-loop NMHV hexagon and fixing previously undermined coefficients in a recent 
Ansatz for the three-loop MHV hexagon. In addition, we consider amplitudes restricted to a 
two-dimensional subspace of Minkowski space, and obtain some equations which involve only 
that sector. 
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1 Introduction 



M = A supersymmetric Yang-Mills theory (SYM) is believed to be integrable in the planar 
limit [1, 2], cf. [3] for a recent review. This has made it possible to compute quantities of the 
theory, such as the spectrum of anomalous dimensions, at any value of the coupling [4-6]. 
On the other hand, remarkable structures in the S-matrix of the theory have been unraveled 
recently. Among them, a hidden, dual superconformal symmetry has been discovered both 
at strong [7] and weak coupling [8] for the S-matrix, which, together with the ordinary super- 
conformal symmetry, generates an infinite-dimensional symmetry encoding the integrability 
of the theory, the so-called Yangian symmetry^. 

The dual superconformal symmetry can be understood as the symmetry of null polygonal 
Wilson loops in a dual spacetime. The bosonic Wilson loops are dual to maximally- helicity- 
violated (MHV) scattering amplitudes at strong [7] and weak coupling [11-14], and recently 
the duality has been generalized to the case of arbitrary helicity (N'^'MHV) amplitudes and 
supersymmetric Wilson loops [15, 16] (or with a closely related light-cone limit of correlation 
functions [17, 18]). Although generally tree amplitudes are Yangian invariant [19-22], the 
naive Yangian symmetry is broken for loop- level amplitudes /Wilson loops even if we consider 
finite quantities, such as the remainder and ratio functions [19, 23, 24]. In this paper, we will 
argue that the Yangian symmetry can be made exact for all- loop amplitudes/Wilson loops, 
which are in turn completely determined by the all-loop equations derived from the exact 
symmetry. 

In order to discuss regulator-independent relations in a uniform way, it proves convenient 
to introduce the BDS-subtracted S-matrix element Rn,k, 

An,k = X Rn,k (1.1) 

where An,fc stands for the N'^MHV scattering amplitude and ^^^^ for the exponentiated 
Ansatz proposed by Bern, Dixon and Smirnov (BDS) [25], including the MHV tree and 
coupling constant factor. 

The BDS-subtracted S-matrix Rn^k is infrared finite and regulator independent. It is 
invariant under the action of a chiral half of the dual superconformal symmetry as well as 
under dual conformal transformations. In this paper we propose a compact, all-loop equation 
for the action of other dual superconformal generators, denoted as Q, in terms of a one- 
dimensional integral over the collinear limit of a higher-point amplitude: 

Q^Rn,k = a Res,=o / (d'l'^n+i ) [Rn+iMi - Rn,kRn+Xi] + cycHc, (1.2) 
Jt=0 ^ ' °- 

where a, ^ = 1, . . . , 4 are momentum-twistor indices and e, r parametrize ^n+i in the collinear 
limit (r being related to the longitudinal momentum fraction), Eq. (3.1). In this equation, 

^This was proposed originally in [9], and the form of the Yangian algebra used in this paper, which takes 
dual superconformal algebra as the level-zero subalgebra, has been discussed in [10]. 
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Figure 1. All-loop equation for planar M — A S-matrix. 

a = a{g^) is one quarter of the cusp anomalous dimension 

a:=-V,^,^=g^-—g^ + —g^ + ..., (1.3) 

known exactly at all values of the coupling g"^ = ^^^^ [4, 26, 27] . We expect the equation to 
be exact at any value of the coupling, but in this paper we will study it perturbatively with 
respect to a. 

We find Eq. (1.2) natural and pleasing in many respects. First, it relates finite and 
regulator-independent quantities. Integrating out a particle with measure {cP^'^Z)^ is virtually 
the simplest operation one could imagine, which carries the quantum numbers of Q. The one- 
dimensional collinear integral over r reflects the physical intuition that naive Q is violated 
because it causes asymptotic states to radiate collinearly. The presence of two terms on the 
right-hand side has a simple explanation: if the first term is viewed as the effect of Q on the 
amplitude, then the term with -Rj^'^x i is due to the action on the ^^^^ factor in Eq. (1.1). 
The proportionality to Fcusp of the second term is thus easy to understand, it being the 
constant of proportionality in the BDS Ansatz, while the structure itself is rigid: certain 
divergences which would violate conformal invariance cancel between the two terms. The fact 
that only 1 — )• 2 splitting appears to all loop orders, as opposed to 1 — )• 3, 4, . . . seems difficult 
to understand from the scattering amplitude viewpoint, and we can only derive it through 
the duality with Wilson loops. 

The equation holds for generic configurations, that is, it neglects so-called distributional 
terms which are supported on singular configurations. These terms were used in [28] to 
determine tree amplitudes. By stripping off the MHV tree, which give rise to such terms, the 
tree amplitudes have been argued to be uniquely determined by requiring analytic properties 
such as the right collinear behavior, in addition to Yangian invariance [23]. In this spirit, we 
will assume that all the pertinent information is included by imposing in addition to Eq. (1.2) 
the correct collinear limits of BDS-subtracted amplitudes, which play the role of boundary 
conditions to Eq. (1.2). 

Using the discrete parity symmetry of scattering amplitudes, we can derive an equivalent 
equation for the level-one generator 

/^(l)«E> ^ay f°° dT{dr]n+l)A r? ST n 9^n,k\ , ,. .... 

Qa ^n,k = a ^nl™ / Rn+i,k " 2^ C'„,ij— ^ + cyclic, (1.4) 

^ \ l<j<i<n-3 I 
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where Cn,i,j is given in Eq. (3.24). The level zero generators Q and Q together with Q^^^ 
generate the full Yangian algebra. 

It is significant that the right-hand sides of Eq. (1.2) and Eq. (1.4) take the form of 
linear operators acting on the BDS-subtracted S-matrix (viewing R^^i as a collection of 
constants, see (2.4)). This means that the right-hand side could be moved to the left, and 
the resulting operators interpreted as quantum-corrected Q and Q^^^ which annihilate the 
S-matrix. In other words, these equations are not "anomaly equations" — their content is 
precisely that all Yangian anomalies can be removed by a simple redefinition of the generators. 
Symmetry generators which receive quantum corrections but nonetheless admit simple closed 
form expressions are not uncommon in integrable systems, see for instance [29]. This being 
said, we will continue to write these equations in the form of a naive (or bare) generator on 
the left, with the correction on the right-hand side, as this will prove most useful for our 
applications. 

The form of (1.2) is very similar to that obtained by Sever and Vieira in the context of a 
proposed CSW-regularization of amplitudes [30]. The essential new features are the focus on 
the finite quantity R, which gives rise to the differenced form, and the advantage of working 
with integrated amplitude, which results in an all- loop relation with an overall proportionality 
to Fcusp- Our formulas also reproduce the one- loop results of [31]. 

Our derivation of Eq. (1.2) will be based on the Operator Product Expansion (OPE) 
for null polygonal Wilson loops [32]. It will be supported by an explicit computation of the 
fermion dispersion relation to ©(F^ygp), finding agreement with [33]. We will also present 
strong explicit evidence for its all-loop validity, through two- and three-loop computations. 
In particular, we have reproduced the two-loop MHV [34] and NMHV hexagon [35], and 
obtained new results for NMHV heptagon and three-loop MHV hexagon, where we have fixed 
the two undetermined coefficients in the recent Ansatz for its symbol [36]. 

The paper is organized as follows. We start in section 2 with a short review on momentum 
twistors, BDS Ansatz and Yangian symmetry. In section 3, we explain how to use (1.2) to 
compute Q, especially at one loop, how are MHV and NMHV amplitudes uniquely determined 
by (1.2), and how to derive (1.4) from it. We then employ the method to reproduce results 
for two- loop MHV amplitudes in section 4. We continue in section 5 to derive results for two- 
loop NMHV hexagon and heptagon, as well as three-loop MHV hexagon. In section 6, we 
outline a derivation of the equation. We also discuss two-dimensional kinematics in section 7. 
We finish with some conclusions and appendices containing techniques and details of some 
computations. 

2 Momentum twistors, BDS Ansatz and Yangian symmetry 

Since our discussion will center on the dual superconformal symmetry, it is advantageous to 
use momentum-twistor variables introduced by Hodges [37] , which manifest the symmetry at 
least for tree amplitudes. The Wilson loop dual to n-point amplitude is formulated along a 
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n-sided null polygon in a chiral superspace with coordinates {x,6), 



xf-xn = AfAf, eT^-e:^, = Xtr^f, (2.1) 

for i = 1, . . . , n, where a, a are SU(2) indices of spinors Aj and their conjugates Aj encoding 
the null momenta of n particles, and A is the SU(4) index of Grassmann variables r]i describing 
their helicity states. The momentum (super-) twistors are defined as 

= {Zt,xt) ■■= {X2,xf\^a,e^^X^a). (2.2) 

We further define the totally antisymmetric contraction of four twistors, {ijkl) := eabcd^fZjZ^Zf, 
and the basic R-invariant of five super- twistors, 

[ijklm]— 5^\^{{{ijklm))) ^2^^ 

{ijkl){jklm) (klmi) (Imij) (mijk) ' 

where the argument of Grassmann delta function is {{ijklm))^ := xf{jklm)+ cyclic. NMHV 
tree (divided by MHV tree), appearing in (1.2), is simply given by a sum of these R-invariants 

l<i<j<n 

At loop level, the symmetry of amplitudes is broken by infrared divergences, which need 
to be regulated and subtracted for exact symmetry. Based on the known infrared behavior and 
the ABDK iterative relation [38], EDS have proposed an exponentiated Ansatz for all-loop 
MHV amplitudes in D = 4 — 2e dimensions [25], 



/iBDS °° 



/I tree 
^n,MHV 



E^'' (/('H^)MW(^e) + cW + i^W(e) 



.£=1 



(2.5) 



where g'^ := 2g'^(47re '^Y has been used as the parameter of loop expansion, f^^\e) = \T^cvls-p + 



0{e), C^^^ are independent of kinematics or n (non- vanishing for £ > 1), and En vanish as 
-)• 0; by stripping off th( 

■■= ^iMH'v/<MHV is given by 



e ^ 0; by stripping off the MHV tree ^''mhv = ^"^^(i^2K23)°.'.'(?i)^''''^ ' one-loop amphtude 



l-loop 



n 



n / 2 \ ~ 

(2.6) 



where Fn^""^{e) is a sum of finite parts of the so-called two-mass easy box functions [25]. The 
BDS Ansatz is believed to be exact for n = 4,5, in which case R^fi = R^^ = R^^i/R^^f = 
1, and R behaves simply under collinear limits, both k -preserving Rn,k — ^ -^n— i,fc s-nd. k- 
decreasing j [1,^2 n-in 1 2] ~^ Rn-i.k-i) as Zn Zn-i (using a parametrization such as 
Eq. (3.1)). 
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As a consequence of Poincare supersymmetry of scattering amplitudes, the BDS sub- 
tracted amplitude is invariant under a chiral half of dual superconformal symmetry and the 
R-symmetry, and it is also believed to be invariant under all bosonic symmetries, including 
dual conformal symmetry, 

i=l ^^i i=l ^^i i=l ^^i 

Naively the BDS-subtracted amplitude is not annihilated by generators in the other chiral 
half, 

Q^ = {StQi):=Y.^fg^, (2.8) 

i=l * 

but as we can see from (1.2), the symmetry is restored by a quantum-corrected Q. Note the 
correction is manifestly Q-invariant, thus the Q invariance and (1.2) imply the invariance un- 
der = ^{Q^, Q^}. For Yan gian symmetry, one needs at least one additional level-one gen- 
erator, e.g. Q^''" which contains the ordinary superconformal generator := Y17=i gx~%^' 

Note that, although second order in derivatives, this is only first order in bosonic derivatives. 
Since is the parity conjugate of s^, which is part of Q, we will derive the Eq.(1.4) from 
the conjugate of Eq.(1.2). 

3 The Q equation 

In this section we elaborate on the evaluation of Eq. (1.2). It involves adding a particle in a 
collinear limit. In the case of edge n, we parameterize its (super-)twistor as 

Zn+i = Zn- eZn-i + CerZi + C'e^Z2, (3.1) 

with C = ^"yti23)^^ ^^"-^ ^' ~ (ra-2n-i2i) " collinear limit is e — )• and, physically, r is 
related to the momentum fraction shared by particle n+1 in that limit. The most general 
collinear limit would require three parameters, and the third one, which we will not need, 
could be obtained by replacing C' with an order one function. The signs and normalization 
have been chosen such that, for real e > and r > 0, Euclidean n-gons (configurations with 
positive cross-ratios) are approached by Euclidean (n+l)-gons. 

The basic operation reSe=o / {d^^'^Zn+ija can be evaluated as follows. In our parametriza- 
tion, the bosonic part of the measure is {(fZn+i)a '■= {Zn+idZn+idZn+i)a = C{n—lnl)aededT, 
where only the dominant part at e — )• was kept. Thus 

rT=oo r- POO 

res,=o / (d^^'^Zn+i)^ = C(n-lnl), res,=o / ^de / drid'^Xn+if- (3.2) 
Jt=0 J Jo 
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The notation res£=o means to extract the coefficient of ^ ™ the e — )• hmit. It is not trivial 
that this exists, but we will see that there are never singularities stronger than ^ log^""*^ e in 
Q of the £-loop amplitude, and that the logarithms always go away after r integration. 

3.1 R- invariants 

It is useful to illustrate the procedure on the simplest non-trivial object, NMHV R-invariants. 
If the R-invariant does not involve Zn+i, the Grassmann integral will produce zero. Further- 
more, even if Zn+i appears, a pole 1/e will be absent unless Z„ is also present. Thus the only 
R-invariants which give non-trivial results contain both and Zn+i- 

Consider the invariant [ij fcnn+1] for z, j, k all distinct from n—1 and 1. After doing the 
Xn+i integration one gets 

/(<!*2„„)^[uin»+ll = C(„-l„l)<. / ^mknn^))Hiik^^^ ,3 3, 

J J {ijnn+l){jknn+\){kinn+l){ijkn+l) 

and plugging in the parametrization Eq. (3.1) and keeping the dominant term as e — )• gives 

de f°° r {{ijknB))^{ijkn) 



pT=00 poo 

reSe=o / {d^^^Zn+i)a['ij knn+l] = C{n-lnl)aVeSe=o— / dr 
Jt=o £ Jo 



{ij nB) (jknB) (kinB) 

-C(n lnl)aj^ '^\^jnB){jknB){kinBy ^"^'^^ 

where Zb ■= Zn~i — CtZi. We could perform the r integral here, but it is advantageous 
not to do so and keep the r-integrand untouched at this stage. This is because in later 
applications we will need this integral with additional dependence on r inserted. However, 
we can simplify it a bit. It has three poles hence two linearly independent residues. Define 
the bitwistor X = X{t) := n A B. Then the residue at (ijX) = gives 

{{ijkn[n-l))^{l]ijn){ijkn) {{n-lnlij))^ 

[n—lnl)a-, r-r-r F ^T~^ w"; r \T~i T — ('^"l^lja^ T- W-^) 

(ijnl) {jkn[n—l) {l\ijn) {kin[n—l){l\ijn) {n—lnli) {n—lnlj) 

This can be rewritten using the nice identity 

n-lnl ^\ =Q^log)^, 

where (n) := (n— Inl). Similarly the residue at (ikX) = is Qlog and the residue at 
{jkX) = is given by minus the sum of the two. Adding these contributions we have, 

res,=o ^^d'\'Zr.^^[^Jknn+l] = T f dlog log M + ^log log . 

Jt=o Jo V {m) {Xik) {m) J 

(3.6) 

This is valid at the level of the r-integrand, for i, j, k ^ n—1, 1. Other R-invariants are com- 
puted similarly, and we complete this subsection by giving the result of the reSe=o Ir=o° ^^'^^n+i 
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operation on these R-invariants, with i,j ^ n—1, 1, 

[ijn—lnn+1] — )• J dlogj^^ — ^ — ^Qlog 



{Xn—2n—l) {ni) 



[ijnn+ll] ^ J ^ log 1^121^ log 



(ni) ' 



^ f dlog^^^-^ Qlog^. (3.7) 

All other R-invariants giving zero. Note these expression all hold at the level of the r- 
integrand. 

3.2 Q of one-loop amplitudes 

Armed with just this result, we are ready to evaluate the Q of any one-loop amplitude. The 
right-hand side of Eq. (1.2) reads 

pT = 00 

QRl;^r = reSe=0 / d'l=^2n+l {RTlMl - RnfRnlXl) + Cyclic. (3.8) 

Jt=0 

Using the (P)BCFW formula for removing Zn+i (associated to the shift Zn+i — >■ Zn+i + 
zZi) [39], the parenthesis can be rewritten as 

n-2 / k \ 

^[nn+llii+1] j;i?*;>(rrri,l,...,i,/OC;Ufc-fc'(^-^+l'---'^)-<fc (3-9) 
i=2 \fc'=0 / 

up to 2^„-|-i-independent terms which do not contribute to the integral. The dependence 
on Zn+i is in the R-invariant and in the shifted twistors, n+1 := (nn+1) n (lii+l), /j := 
(M+l)n(nn+ll). However, the parenthesis has a smooth collinear limit since no term depends 
simultaneously on both Z„ and Zn+i- Thus we can merely replace Zn+i by Zn and by its 
limit li = (n— Inl) H (ii+l) (supersymmetrically, the fermions of are taken from Xi and 
Xi+i)- The Zn+i dependence is limited to the R-invariant and using Eq. (3.7) the integral 
gives 

QRn'r = E rflog ^^^Q log (parenthesis in Eq. (3.9)) + cyclic. (3.10) 

Jt=q ^ (-'^12) (m+1) 

In the parenthesis, n+1 n,ii ^ li, and nothing depends on r. 

We must show that this integral is convergent at its endpoints. Near r = 0, there is a 
pole due to dlog{Xn—2n—l) in the term i = n— 2. However, the two terms in the parenthesis 
cancel in this case, so there is no problem. There is also a pole at r = oo, due to dlog(X12) 
present in every term. It is nontrivial to see that it cancels out in the sum, but can be proved 
as follows. 

Instead of using the (P)BCFW formula associated with shifting Zn+i, we could have 
used the BCFW formula associated with the shift Zn — )• Zn + zZn~i- Then we would have 
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obtained the same formula but with the R-invariants replaced with [ii+ln—lnn+1], and 
so dlog (^xn-2n-i) integrand, but the parenthesis in the e — )• limit unchanged. This 

form would make convergence at r = oo manifest but not at r = 0. We conclude that overall 
convergence follows beautifully from the equality of the BCFW and (P)BCFW representations 
of tree amplitudes. 

Given these cancelations, we can integrate Eq. (3.10) termwise by dropping the dlog(X12) 
factor and the i = n—2 term, obtaining simply 



n— 3 

/=>d1-1oop >r^i (nlu+1) ^ (ni) 




(n, 1, . . . , i, Ii)R'"'%Ii, . . . , n) - <T + cyclic, (3.11) 



which agrees with the formula of [31] (there the product of tree amplitudes is interpreted 
in terms of unitarity cuts). The subtraction of R^'j^ arises because we are considering the 
BDS-subtracted amplitude, which is the one-loop N'^MHV ratio function in this case. 

3.3 Uniqueness of Q solutions at MHV and NMHV 

The Q equation is especially interesting because as we will see now, it fixes uniquely MHV 
and NMHV amplitudes (assuming that the right-hand side is known). This is not too difficult 
to see for MHV amplitudes using the momentum-twistor form of Q, (2.8). Indeed, taking 
derivatives of the equation Qf{Z) = for any function of bosonic Z's, f{Z), we have, 

^,Qlm = ^ ^JiZ) = 0. (3.12) 

This equation, for all particle labels i and twistor indices a = 1 ... 4, implies that a bosonic 
function annihilated by Q is a constant. Thus the ambiguity of the Q equation is at most 
a constant, which can be fixed using the properties of the BDS-subtracted amplitudes in 
collinear limits. 

For NMHV amplitudes, we have to work harder to restrict the kernel of Q. A simple 
example which illustrates this at 5-points is [12345] log |j|||y- This has vanishing Q because 

rio...i/M„ Jl234) ^^^^^^^ (123)((12345)) 

(1235) = (1234)(1235) ^^'^^^ 

contains ((12345)) both explicitly and from the Grassmann delta function 5^l^(((12345))) 
in the R- invariant, hence vanishes. On the other hand, this expression is not acceptable 
because the argument of the logarithm is not conformal invariant (synonymous with little 
group invariance in what follows): it has non- vanishing weight with respect to 4 and 5. So 
it does not correspond to any real ambiguity. This turns out to be general: any NMHV 
expression with neutral little group and annihilated by both Q, Q, is a sum of R-invariants 
with constant coefficients. 
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To prove this, we first note that by Q invariance alone, any NMHV expression can be 
written as 

F= Yl [ljklm]Fj^k,i,m{Z) (3.14) 

2<j<k<l<in<n 

where the ^ ^ ^ [Ij/c/mJ's form a basis for all independent NMHV R-invariants at n- 

point [40]. Each Fj^k,l,miZ) is a conformal invariant function of the bosonic Z's. 

To show that the Fj^k,l,miZ) must be constant, we pick {l,j,k,l,m} and extract 
a specific component, xjx]x'kXfxfm of Z'-Q\F. The only way x] can arise is either from 
QFj,k,l,m or from a R-invariant, but since Zj-r^F = 0, it can only arise from a R-invariant. 
Since only [1 j kl m] contains the prescribed components, we deduce that 

^i^^.>,^- = 0- (3-15) 

Repeating this with permutations of j,k,l,m shows that Fj^^^i^^ is independent of Zj, and 
repeating for other z's shows that Fj^^^i .^ depends only on twistors 1, j. A;, I, m. But since there 
are no nontrivial little-group invariant functions of five twistors, Fj^k,l,m must be a constant, 
QED. 

All remaining constant ambiguities can be fixed by collinear limits. As mentioned, there 
are both fc-preserving and fc-decreasing collinear limits. It turns out that just four of the k- 
preserving limits suffice for any n. For instance, working in the same basis, the /c-preserving 
collinear limit Zi — )• Z2 will fix all constants except those multiplying invariants of the form 
[12...]. Taking the limit Z2 — )• -Z3 will then fix the coefficient of all but those beginning with 
[12 3...], and so on. 

These results open up the possibility of using the Q equation to compute nontrivial 
MHV and NMHV amplitudes. Given one-loop N'^MHV amplitudes as the seed for recursion, 
this will restrict the applications in this paper to two-loop MHV and NMHV and three-loop 
MHV. To go beyond NMHV it becomes necessary to use both the Q and Q^^^ equations^, 
or, equivalently, the Q equation and parity. Uniqueness then follows from a theorem proved 
in [41, 42]: all Yangian invariants are combinations of compact contour integrals inside the 
Grassmannian G{k,n). We conclude that any N'^MHV expression annihilated by (naive) 
Q, Q, Q^^^ can only be a combinations of such invariants, multiplied by c-numbers, which we 
expect to be determined by collinear limits. 

3.4 The one-loop NMHV hexagon 

In the case n = 6, equation (3.11) evaluates more or less directly to 

Q<1°°P = ((5) + (3)) log U3Q log 1^ + (1) log nsQ log |^ + cyclic, (3.16) 



simple counter-example to Q uniqueness is the invariant ^ ' ^|i234)-^(6i2"3)^'^'"^^ which arises in the 6-point 
N^MHV tree amplitude and depends on six twistors. Any conformal invariant cross-ratio of the six twistors 
multiplying it will be Q-invariant. 
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where (1) is the R-invariant [23456], (i) is obtained by a cychc shift, and we have used that 
R^i^ = (1) + (3) + (5) = (2) + (4) + (6). The appearance of log U3 is easy to understand from 
the two terms in Eq. (3.11), because they correspond to two poles of the r-integral, and so 
what multiplies the log has to be equal and opposite. We use the following cross-ratios 

_ (1234)(4561) _ (2345)(5612) _ (3456)(6123) 

~ (1245)(3461)' ~ (2356)(1245)' ~ (3461)(2356) ' ^ ^ 

We have just shown that the information in Eq. (3.16) should suffice to determine -Re.i- 
Let us see how this works. The crucial step is to bring the right-hand side to a form where 
the argument of each Q is the logarithm of a conformal invariant cross-ratio; this form will be 
unique. This can be achieved by adding suitable combinations of zero in the form of equation 
(3.13), for which there is a systematic procedure. 

The following procedure reduces this to a simple linear algebra problem. The first step is 
to remove the ambiguities in writing the R-invariants, by using the identity (1) — (2) + (3) — 
(4) -1- (5) — (6) = to remove (6). We can then use four distinct nontrivial representations of 
zero, [(1) - (2) + (3) - (4) + (5)] times Q^^^, Q^^^, O^fff, to remove e.g. the 

little group weight with respect to i of the coefficient of (i), for i = 1, . . . 4. 

Actually, there is a final constraint: it is not trivial the little group weight with respect 
to 5 of the coefficient of (5) is also removed; but this is is the case. Then the coefficient of 
(i) has correct little group weight with respect to i for i = 1, . . . ,5. The little-group weights 
with respect to other variables can then be removed using equation (3.13) with R-invariants 
(1),...,(5). 

This procedure is simple to follow but not particularly illuminating, so we spare the 
reader the details, recording only the final result: 

<3<1°°" = (<f Olog - ((1) + (4))Qlogn2 - ((2) + (5))Qlogni) logug + cyclic. 

(3.18) 

This equation is equal to Eq. (3.16), but now the Q acts on conformal invariants. The upshot 
is that in this form we are allowed to directly integrate Q: 

RIT" = <"(log^2 log^/3 + Li2(l - u^)) - ((1) + (4)) logns logus + cyclic + C7, (3.19) 

where C is an undetermined combination of R-invariants with c- number coefficients. 

To fix C, we can consider collinear limits. For instance, the ratio function should vanish 
in the /c-preserving limit Zq ^ Z^, corresponding to ui — )• and U2 — )• 1 — U3. This limit 
probes the coefficient of (5) plus the coefficient of (6). In this limit, what we have in Eq. (3.19) 

2 

goes to ^[12345]. All other /c-preserving limits go to the same number, allowing us to fix 

C = -y<f- (3-20) 

This is the correct ratio function! 
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3.5 Derivation of equation (1.4) from equation (1.2) 

This subsection lies a bit outside the main scope of this paper. As noted in the Introduction, 
the equation for Q together with parity symmetry of scattering amplitudes implies an equation 
for 

To derive it, the first step is to express Eq. (1.2) in the language of scattering amplitudes. 
We only need to do this for the two components of Q which coincide with the ordinary 
superconformal generators [9] . Technically, we really only need to do this for the first term 
in the parenthesis of Eq. (1.2), and we can drop the explicit dependence on e, reinstating it 
at the end. We find, after reinstating the MHV prefactor and changing variable Ct — )• 

s;^^n,fc = A„a lim / dx(d^x)'^>in+i,fc+i(---,{A„,xA„,XT/„+x},{A„,,(l-x)A„,(l-a;)r/„-x})+... . 

The variable x is the usual longitudinal momentum fraction. With the help of the BCFW 
computer package for the evaluation of tree amplitudes [43] , we have verified that the integral 
gives the correct result acting on the NMHV 5,6,7 point tree amplitudes. The upshot is that 
in this form it is possible to immediately write down the parity-conjugate equation:'^ 

s'X^n^k = K. / ^^77-^ A+l,fc(- ■ ■ , {A„, XA„, Xfjn+X], {An, (l-x)An, {1-X)fln-X]) + - ■ ■ , 

where the denominator l/x(l — x) comes from a little group transformation needed after 
interchanging A and A. The final step is to convert this equation back to momentum twistors: 

Q'^I^" An,k = Z^lim f —{dxn+i)ARn+i,k{l,-.-,n+l) + ... (3.21) 

where we have put back the e dependence, -Z„_|_i(e, r) being again given by Eq. (3.1). Strictly 
speaking, gives two out of the four twistor components of Q^^". The remaining two 
components come for free, because the level-zero conformal symmetry of Wilson loops is 
unbroken acting on BDS-subtracted amplitudes. 

This takes care of the first term in the parenthesis of Eq. (1.2). To deal with the second 
term we need the explicit form of acting Q on NMHV tree. 



n—3 _ 

reSe=oa Zn+iK^+n - > . d log ^ — . (3.22) 

' ' ^— ' I n 0^ I n \ \ In •n-\- \ n — An — l\ 



((n^^+l)) (nn+lii+l) 

r"log r. 

{ni){ni+l) (nn+ln— 2n— 1) 



^The easiest way to derive this equation is to consider the case where particle ti is a positive-helicity gluon 
in the equation. Then the x integral gives eight terms on the right hand side, involving a minus-helicity 
fermion, or a scalar plus a plus-helicity fermion, with various R-symmetry assignments. Parity dictates that 
when n is a negative-helicity gluon, s'X should produce the eight parity conjugate terms. The correctness of 
the other cases follows by supersymmetry. 
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We take its parity conjugate using Xi = Yl]=n-2i'^j)^j then r]j — )• In terms of 
momentum twisters, the end result is 

QA^'^Rn^k = a K^'^J^ [ — i{dXn+i)ARn+i,k - ^ri,i,j^^^^ ]+ cycHc (3.23) 

-^0 \ l<j<i<n-3 / 



where _ 

{ni){ji+l) - {ni+l){ji) d (nn+lu+1) , . 

(-^n,i,j{r) = / A/ ^TT\ 7 n TT, (3.24) 

(m)(m+l) ar (nn+ln— zn— 1) 

which is the formula recorded in the Introduction. 

Because this is a consequence of unbroken parity symmetry and the equation for Q, this 
does not require separate verification. To ascertain that Eq. (3.23) contains no mistake, we 
have tested it on known expressions for 1-loop 6,7,8-point NMHV amplitudes. 

4 Two-loop MHV amplitudes 

Armed with just the Q equation (1.2), expressions for one-loop NMHV ratio functions, and 
the (P^^Z integral of R-invariants Eq. (3.7), we are now ready to analyze two-loop MHV 
amplitudes. 

4.1 The square and the pentagon 

In the cases n = 4 and n = 5, it is well-known that logi?4^5 = [25]. Let us begin by 
reproducing this simple result starting from equation (1.2). For n = 4, this is essentially 
trivial for all loops provided log i?5 = at one lower loop order, so the equation reads 

QR4 = 0. (4.1) 

This implies that -R4 is a constant, which must be trivial by the boundary condition. 

In the case n = 5, starting at two loops, the right-hand side is not so trivially zero. Rather, 
it involves the collinear limit of the one-loop six-point NMHV amplitude given in Eq. (3.19). 
Letting 

res,=o <f\'Z,R]:'^"^ = Q log x /, (4.2) 

Jt=o (4513) 

we get that the R-invariants contribute to / as follows, 

(l)^dlog^, (2)^dlogT, (4)^-(ilog(T + l), (3),(5),(6) ^0. (4.3) 

In the collinear limit ui — )• e^, ii2 — )• jj^^ — ^ allowing us to write I = Ii + I2 with 
h= d (log \ ] log(l + r) log(l + 1/r), 



/■oo 

h = log X / d (log(l + r) log(l + l/r)) . 
^0 



(4.4) 
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We find that Ii = I2 = 0, confirming tliat QR^ = as expected. This is the first nontrivial 
hint that the equation is working beyond one-loop. 

The reader might worry about the divergent prefactor loge^ in front of I2. Shouldn't 
the e — )• limit entering our basic equation be well-defined? The answer is that the order of 
operations is important. The limit e — >• will always be well-defined provided the integration 
over r is carried out first. If one were to take instead e — >■ with fixed r, one would find a 
divergence. This divergence has a simple explanation and is actually predicted by the Wilson 
loop OPE [32]. We will return to it in subsection 6.3 where we confirm the quantitative 
prediction for it, and also give the general argument for its cancelation after r-integration. 

4.2 The hexagon 

For n = 6, we need the one-loop seven-point NMHV amplitude, which can be put in a compact 
form [44], 

= [1, {2, 3}, {4, 5, 6}]/i + [1, {2, 3}, {4, 5, 6, 7}]/; + cyclic, (4.5) 

where 

[z,{m,...,j},{A.,...,/}]= E [^"^'^]' (4-6) 

J={i+l,i+2} K={k,k+l} 



h = Li2(l - v^vz) + Li2(l - vi) + Li2(l - v^vq) + Li2(l - 

- Li2(l - viVi) - Li2(l - vq) - Li2(l - v^) - Li2(l - v^vi) + log7;7logf2, 
I[ = Li2(l - v-j) + Li2(l - vq) + Li2(l - V3) + Li2(l - v^vi) 

- Li2(l) - Li2(l - v^v■i) - Li2(l - V2.vq) + log^uylogfe, (4.7) 

and are obtained by cyclic shifts for i = 2,..., 7. Here we need to define cross-ratios 
beyond six points 

_ {ii+ljj+l){kk+lll+l) 

and at seven points, a basis of cross-ratios can be chosen as Vi := iii+i,i+3.i+4,j for i = 1, 7. 
In the collinear limit — )• Zq, — )• 0, ^3 — (1 — V2)/{1 — V2'V6) and ^ (1 — vg)/{1 — V2Ve), 
thus the result depends on vi,V2,vq,V7. 

The R-invariants appearing are not independent, and it is convenient to choose those 
containing the label 2 as a basis. Upon doing the integral over cP^^Zj, only [12367], [12467], 
[23467], [23567] and [24567] contribute ([12567] does not contribute, because its coefficient 
has to vanish due to the k-decreasing collinear limit constraint), which produce, proceeding 
as in the five-point example, 

QRl-^r = ihA + h,2)Q log 1^ + (l2,i + h2)Q log 1^ + cyclic, (4.9) 
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where 



ceo 



'1,1 



'1.2 



(ilogi 



l0g^^^l0g^^ + Li2(l-n3)-Li2(i-^)) \ 



+dlog(r + l) (^logui(r+l)log^+Li2(l 



U3) 
Us) 



Li2(^ 



■+"3 ' 
-■"3) 
+U3 



V 

loge^ X 



+Li2(l - ui) + Li2(l - U2) + logui logU2 - ^ 



n3(T + l) 
r + Us 



log( 



T + U3 



+ log(T + 1) log 



T + ^3 

r+1 



and 



/ dloe 



2,1 



^2,2 



JO 



T+H3 



log ^^^^ log - log(r + 1) log 

+Li2(l - M2) + Li2(l - U3) - Li2(l - 



r+l 



(4.10) 



= loge X 



+dlog^(Li2(l 

T + U3 



+dlog{T + us) log ^ log ^ + Li2(l - Ui) - Li2(l 



^)-Li2 



"1 

T+U3 



■ l-"3 
^ T + 1 



T+«3 



+ Li2(l-^)+log^log^ 



d log 



r 



log 



^^3 



T + U3 



(4.11) 



The non-spacetime ratio v = |5623)(6i24) needed to produce a parity-odd part. 

We emphasize that this comes directly out of the collinear limit of the heptagon. No 
manhandling has been applied, nor would have been necessary. We have, in the interest of this 
presentation, used standard dilogarithm identities to simplify the expression and hopefully 
make it more human-readable, but we have not used integration by part nor any manipulation 
which would affect the numerical value of the r-integrand. 

The divergent terms cancel upon integration: /i^2 = -^2,2 = 0, just as in the pentagon 
example. This cancelation is of paramount importance to our approach, and after it is 
effected, we are left with two finite and manifestly conformal-invariant integrals Ii^i and /2,i- 
The mechanism for this cancelation is general and detailed in subsection 6.3. 

The integrals produce trilogarithms. Computing them is not entirely trivial (for instance, 
Mathematica would not do them automatically), but obtaining their symbols is, following, 
for instance, the method of Appendix A. From the symbol it is not too difficult to obtain 
actual functions, and then fix beyond-the-symbol ambiguities using the differential computed 
in Appendix A. The resulting functions are quite simple 



I2A 



^ logins + logui logU2 + ^Li2(l - Uj)^ log'U3 - 2Li3(l - 1/^3), 

- + Li3(l - I/U2) + Li3(l - l/us) - Li3(l - 1/ui) + ^ log3 ^ 
2 12 ui 

, 1 1 U2U3 sr^j. . , ^ 
+ -log 2^Li2(l - 1/Ui), 

^ ^1 i=i 



(4.12) 
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where Iq^ is the six-dimensional massless hexagon integral [45, 46], reproduced in Appendix B 
alongside the definitions of and to be used shortly. 

To complete the computation of QRofl = dRe^, we go back to Eq. (4.9) and add the other 
edges contribution via cyclic symmetry. For future reference, we record the simple result: 

di?2J°P = Il^dlog — + f /i,idlog + two cyclic) . (4.13) 

This agrees precisely with the differential of Goncharov, Spradlin, Vergu and Volovich's for- 
mula [34], derived from the results in [47], 

3 / 3 \ ^ 

i=l ^ ^ \i=l / 

(4.14) 

Of course, in practice the step from Eq. (4.13) to Eq. (4.14) can be a very difficult one, 
and we do not wish to imply otherwise; we have simply gone the other way, taking the 
derivative of Eq. (4.14). Still, it is impressive how close to Eq. (4.14) the present formalism 
lands us, namely, on Eq. (4.13). Important qualitative features of the result, such as its 
finiteness, transcendental degree and conformal invariance, were manifest at every stage of 
the computation. 

4.3 The differential of the n-gon 

To obtain results for n > 6 two-loop MHV amplitudes, we need the (n+1) > 7-point one- 
loop NMHV amplitudes. Since there are no qualitative differences between (n+1) > 7-point 
amplitudes and the seven-point one, the computation is similar in every respect. We have 
verified that the divergent terms integrate to zero for generic n, leaving a set of finite and 
manifestly conformal integrals, which are too lengthy to record here. However, we have 
explicitly obtained these integrals for n = 6, 7, 8, 9 (which is generic) using the present method, 
and we can compare this result with that given in [24] (specifically, equations (4.21) and 
(4.28) there). We find perfect agreement: numerically both one-dimensional integrals give 
the same to 30-digits precision on a few randomly generated Euclidean kinematic points, and 
symbolically, they give the same symbol. We recall that these integrals give degree-three 
transcendental functions characterizing the full differential of the amplitudes. 

5 Two-loop NMHV and three-loop MHV amplitudes 
5.1 The two- loop NMHV hexagon 

Because one-loop N^MHV amplitudes are known, there is no reason to stop at MHV level. 
(For one-loop amplitudes, we have used expressions based on the box-expansion and gen- 
eralized unitarity expressed in momentum twistor space [22, 48-51].) The first step in our 
procedure to compute the NMHV hexagon is to take the collinear limit of the one-loop seven- 
point N^MHV amplitude and extract the de/e term from the d'^^^Z integration. Just as in the 
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previous cases, one obtains a one-dimensional integral over the variable r, and after verifying 
that terms proportional to log e integrate to zero, one is left with a manifestly finite and con- 
formal integral over polylogarithms of degree two. The integrals are not significantly more 
difficult than those appearing in the MHV case, and can be done similarly; we only record 
the result: 



= (6)0log||/i + ((1) - (2) + (4) - (5))Qlog||/2 + ((2) - (4))Qlog||, 



+ ((6)Q log II + ((5) - (4))Q log /4 + ((2) + (4))Q log ||/5 

+(5)Qlog^/6 + (3)QlogM/,, (5.1) 
(63) (63) 

where /i, . . . , /y are degree-three transcendental functions reproduced in Appendix B. 

The fact the result could be expanded over a basis of 7 linearly independent rational 
prefactors (of the form (R- invariant) xQ log l^jy), times integrals with unit residues, follows 
from a general Grassmannian analysis. The key fact is that these prefactors all originate 
from seven- point N^MHV leading singularities, which are combinations of the 15 independent 
residues in the G{2,7) momentum twistor Grassmannian [21, 22]. In fact, we found that 
coming up with the full list of the 7 prefactors was the most nontrivial part in our derivation 
of the above equation. After this was known, the reSe=o part of the (P^^Z integration step 
could be easily automated on a computer.^ There remained only the r integration, which 
could be done automatically at the level of the symbol and with a bit of human input for the 
function. 

After obtaining this equation, we are (already) essentially done. To complete this compu- 
tation, we need to use cyclic symmetry to obtain the contribution of other edges, and plug the 
result into the exact same linear algebra problem as encountered in the one-loop example in 
subsection 3.4. Namely, starting from the 42 rational prefactors obtained from symmetrizing 
the above 7 over the 6 edges, we need to add "zero" in the form of equation (3.13) to make the 
argument of all Q's become cross-ratios. Just like at one-loop, we found exactly one potential 
obstruction, which vanished for the above /j's, leaving 41 truly independent functions. We 
expect this counting to be the same at all higher loop orders. Expressing the amplitude in 
the form [52] 

Rl'!r = \ ([(1) + (4)]^3 + [(2) + (5)]l/i + [(3) + {^2 

+ [(1) - {A)]V, + [(5) - (2)]Fi + [(3) - (6)]V2) , (5.2) 



then the solution yields the differentials of each F's and V^s. The resulting formulas are 
reported in Appendix B, together with the definition of y variables. 



*We used a semi-numerical procedure, in which we evaluated numerically the d?^^Z integral of the N^MHV 
residues for a set of random integer-valued momentum twistors. We then used the analytic knowledge that 
the result should be an integral linear combination of 7 basic objects to promote the numerical result to an 
analytic one. Although semi-numerical, this procedure has no error bars and is rigorously exact. 
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From these differentials we can already check that the symbols of V and V agree with 
those obtained recently by Dixon, Drummond and Henn [35], attached with their arXiv 
submission; they do. This is one first nontrivial check. But we are also interested in beyond- 
the-symbol information. We could in principle compute the differential of the results in [35] 
and compare with Appendix B, but we have contented ourselves with a numerical comparison. 

To obtain numerical results we first need the value of the functions Vi, V2, V3 and Vi, V2, V3 
at at least one point. Using the fact that V1+V2 and V3 vanish in the ui — )• collinear limit, for 
instance, we could in principle evaluate these combinations at any point by integrating along 
a path connecting to this limit (choosing a path which remains in Euclidean kinematics). We 
would then use other paths to compute the other cyclically related combinations. However, 
we did not find this approach particularly convenient in practice. A more fruitful strategy is 
to first derive the amplitude at some other point away from a collinear limit. In fact, in the 
special case ui = U2 = U3 = u, it turns out that the differential simplifies dramatically 

KD 4 o 47r^ 

dV = -I^^dlogy + {2Lh{l - u) + ALhil - 1/n) + 5 log nLi2(l - u) + - log^ n - — logu)dlogu 

-(6Li3(l - n) + 6Li3(l - 1/n) + 61ognLi2(l - n) + 21og^ n - 27r^ log n)(i log(l - n), 
dV = 0. (5.3) 

where V := Vi = V2 = V3 and V := Vi = V2 = V3 = 0, allowing it to be integrated explicitly 

1 vr^ 1 

V{u, u, n) = -4L+(x+n, x-u) - r-:J^ - + 2{Ui{u) + - log^ u log(l - n)) 

18 9 o 

-6Li4(l - n) - 6Li4(l - 1/n) + 4Li3(l - u) \ogu- 5Li2(l - n)Li2(l - 1/n) 

7 57r^ TT^ 

+ - log^ n - 27r2Li2(l - u) - — log^ u - 2C(3) logn + — . (5.4) 
z4 D ID 

The first three terms are essentially as in — l/3i?g^Q. To fix the constant, we have used 
numerical integration as explained in the previous paragraph, connecting these configurations 
to a collinear limit. We have computed the value of the constant at the three points n = 
1/3,3/4 and 5/6; each point produced the same result. We then recognized this numerical 
result as tt^/IO and confirmed it to 40 digits. Because it is fully manifest from the formulation 
that the constant is a degree four transcendental number with order one rational coefficient, 
it does not seem necessary to supplement the numerics with an analytic computation. 

The upshot of the formula is that it is very easy to deform any kinematical point to 
one on the line ni = n2 = n3 = u. Integrating the differential along such paths, we can 
evaluate efficiently the ratio function at any point. In particular, we have evaluated it on the 
kinematic point in [52]. Defining V' , V by adding the one- loop shift —^Rq[ and multiplying 
by 1/4 to account for expanding in a as opposed to 2*7^, we find on the kinematical point 
(ni,n2,n3) = (^, ff , f) (see Eq. (3.17)) 

V{ = 12.6138748750304719319, V{ = -0.121176561122269858950i (5.5) 

V2 = 11.7057979933899946922, V2 = 0.030638530205807842307^ (5.6) 

= 14.4289552936316184920, = 0.090538030916462016643i. (5.7) 
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This was obtained by integrating along a simple linear path in cross-ratio space to the point 
u = 9/4, but we have also tried a few other points and got the same result. The parity even 
objects Vl agree precisely with the quantities called Vi + in [52] and the parity odd objects 
V- agree within numerical accuracy with those given in [35]. After accounting for the same 
coupling constant shift, Eq. (5.4) can also be compared directly with Eq.(6.30) of [35]; we 
have compared the value at the two points given in Appendix D of [35] and found perfect 
agreement. Given that their symbols match, these numerical tests remove any doubt in our 
mind that the two expressions are equal. 

5.2 The two-loop NMHV heptagon and the three-loop MHV hexagon 

The NMHV heptagon can be attacked in an entirely similar way starting from the collinear 
limit of the known one-loop N^MHV octagon. 

The first step is essentially kinematic and independent of loop order: one has to list 
all (rational) objects which can arise from taking residues of the (P^^Zj integral on octagon 
leading singularities. We found 42 linearly independent ones, all of the form (R-invariant) 
times Q log |?^y , where i and j are momentum twistors or intersections of the momentum 

twistors entering the R-invariants. An example being [23457] Q log ^^^^^jrl^^, but actually 
only three elements of the basis contained intersections. In general at i-\oop we expect to find 
the same 42 structures, each multiplying a pure transcendental integral over degree 2(i — 1) 
functions (with potentially a finite a number of additional ones related to 8-point leading 
singularities not visible at one- loop, which we have not considered). In the case at hands, 
over dilogarithms. Another purely kinematic step is the analog of the linear algebra problem 
encountered previously: out of the 7 x 42 = 294 residues obtained by cyclic symmetry, one has 
to find all combinations which can be written Qlog of (conformal invariant object), possibly 
adding zero in the form of Eq. (3.13). We found 288 combinations, leaving 6 constraints on 
the integrals. These 288 combinations are independent of loop order. 

We have not computed the resulting 42 integrals (each of which, manifestly, would give 
trilogarithms), but we have computed their symbol. This was essentially automatic using the 
method of Appendix A. Plugging the result into the solution of the linear algebra problem 
then gives the symbol of the amplitude. All entries of the symbol are either four-brackets 
or intersections of the type (12(4) n (6)) or (23(745) n (7)). We hope to analyze it further 
elsewhere.^ 

In this paper, our interest in the heptagon stems mostly from its connection with the 
three- loop MHV hexagon via the Q equation. In fact, as already familiar from our analysis of 
the two-loop MHV hexagon, in an appropriate basis out of the 15 independent R-invariants 
at 7-points only five survive d'^^^Z-j integration, namely, [12367], [12467], [23467], [23567] and 

^In its present unprocessed form, the result is too lengthy to be attached with this arXiv submission. It is 
available upon request to the authors. 
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[24567]. If we write cIRq'^q"'^ = dlog |§§x|yA +dlog |||j|y + cyclic, it follows that we can write 
h = [ ( dlog(r + l)gi + dlog ^^^^g^ + dlog ^ ^ g3 ) 

Jo \ r T + U3 J 

h = I [ d\og{T + v)g2 + dlog 55 + dlog ■ — 53 (5.8) 

Jo \ T T + V J 

where v = |5g23)(6i24) • "^^^ functions gi are pure degree- four transcendental functions 
determined by the collinear limit of the heptagon ratio function. On physical grounds (the r 
integrals must converge), we know that gi^2 must vanish at r = 00 and 5(4^5 must vanish at 
r = 0. Thus we can use integration by parts: 

h = [ ( log(r + l)/ii + log ^^^^/i4 + log ^^^^^/i3) +53(0) log n3, 

Jo \ T T + U3 J 

I2 = r flog(r + v)h2 + log — /15 + log ^^h^) + 53(0) log -. (5.9) 

Jo V ^ T + V J U3 

where hi = —^gi are degree 3 functions. We see that only the collinear limit of the differential 
of the heptagon is needed, with the exception of 33 (r = 0), but one could argue that it is 
fixed by cyclic and parity symmetry of dF^^Q°^ . We hope to use these equations in the future 
to study the differential of the three-loop MHV hexagon, beyond the symbol. In any event, 
our result for the symbol of the heptagon already gives the symbols of the gi, which, after 
the nontrivial but entirely automated integration in Eq. (5.8), give the symbols of Ii and I2, 
which in turn give, directly, the symbol of i?g'Q°°''. We now describe this result. 

Recently, an Ansatz was constructed for the three-loop hexagon, based on reasonable 
physical assumptions about entries of its symbol [36] (most significantly, that they should all 
be products of momentum twistor four-brackets), on OPE constraints [32], and on requiring 
that that the last entry of the symbol should involve only brackets of the form {i—lii+\j). 
This Ansatz contained many coefficients but, remarkably, in the end all but two could be 
determined by these authors. Recently Lipatov and collaborators, considering Regge limits 
using new results on the adjoint representation BFKL kernel, confirmed the value of a number 
of these coefficients [53]. 

There are three things we wish to add here. First, that all entries of the symbol should be 
four-brackets is manifest from our approach, since it follows from the symbol of the two-loop 
heptagon involving only momentum twistor intersections (together with the way symbols of 
integrals are built using e.g. the algorithm in Appendix A, and the fact that at six points all 
momentum twistor intersections become reducible to four-brackets). In turn, this property 
of the heptagon was essentially inherited from properties of the one-loop octagon in collinear 
limits. In this way it should also be possible to obtain general information about the symbol 
at £ > 4 loops, although we will not do so here. 

Second, the assumption about the last entry, conjectured in [24], can actually be derived 
from Eq. (1.2) and is therefore now proved. Indeed it follows from writing the NMHV hep- 
tagon in the form (R-invariants) times (pure transcendental functions), and using the general 
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result for (P^^Z on R-invariants, equation (3.7). The upshot is that these two assumptions 
made in [36] follow rigorously from Eq. (1.2), without doing any explicit computation. Note 
that although this form of the heptagon is not strictly proven to all loops, it is widely be- 
lieved that it does hold [21], and assuming this then the statement about the last entry of 
the symbol for MHV amplitudes follows to all loops. 

Third, we have found that our first-principle computation of the three-loop hexagon 
symbol is consistent with the Ansatz of [36], which gives a highly nontrivial check on both 
our approaches. Our new result can be summarized very succinctly: the final two coefficients 
in [36] are ai = — | and a2 = 

6 All-loop validity of the Q equation 

In this section we would like to explain how we believe Eq. (1.2) could be proved, and show 
its consistency at any value of the coupling. 

6.1 Outline of a derivation 

Our proposed derivation of Eq. (1.2) starts from an expression in [16, 24] for the right-hand 
side of Q in terms of insertion of a fermion operator on the edges of the chiral Wilson loop 
(defined on {x,9) space): 

Qi{Wn,k)^9^ fdXaa{{tl^'^ + Fe^ + ...rWn,k). (6.1) 

This can be decomposed into a sum of n terms, one for each edge. Since each edge contribu- 
tion is gauge invariant and meaningful, for the following discussion it will suffice to consider 
the contribution of edge n. For simplicity we will also assume that the (unbroken) Q super- 
symmetry has been used to set fermions Xn-i, Xn and xi to zero. Then ^ = along that 
edge, and the formula reduces to the supersymmetry transformation law of a bosonic Wilson 
line in a suitable normalization. 

The above equation holds for the Poincare supersymmetries of the Wilson loop. For the 
superconformal generators S*""^, extra terms are expected due to the breaking of conformal 
invariance. On the other hand, there is no need to study S explicitly because in the end when 
we obtain equations for Rn, which is known to be conformally invariant. On Rn, the action 
of S is simply related to that of Q. 

In [24], the chiral Wilson loop with fermion insertion was computed in explicit examples 
using conventional Feynman diagram techniques. The new ingredient in the present paper 
is a simple yet powerful fact about the spectrum of excitations of the null Wilson loop: the 
fermion insertion is the unique twist-one excitation with the quantum numbers of Q. 

This is a powerful statement because it means that the right-hand side of Eq. (6.1) isn't 
really a new object. Rather, in the spirit of the Operator Product Expansion (OPE) of null 
polygonal Wilson loops [32], its expectation value can be extracted from any object having 
a nonzero overlap with it in the OPE limit. The rest of this derivation will thus be based 
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Figure 2. Fermion insertion on the Wilson loop versus kink insertion 

on the analysis of [32]. The simplest possible object is the collinear limit of a (n+1)— point 
Wilson loop; this is depicted in Fig. 2. A good strategy to extract the piece with the right 
twist and quantum numbers in this limit is to write down the simplest operation with the 
quantum numbers of Q, namely the (P^^Zn+i operation detailed in section 3. 

To be more precise, the fermion insertion is part of a one-parameter family of insertions 
having bare twist one (at weak coupling), labeled by a position r along the edge. In the 
quantum theory, operators in this family will renormalize among themselves. Thus we expect 
a relation of the form 

/>oo 

limo!2|3z„+i(r,e)(W^„+i,fc(r,e)) = / dr' F{t,t' , e){^P{T')Wn,k) , (6.2) 
-^^o Jo 

where the inserted twistor Zn+i{T, e) is parametrized as in Eq. (3.1), and the right-hand side 
contains the Wilson loop with insertion we are interested in. Now we could instead consider 
the BDS-subtracted Wilson loop, and using the collinear limit properties of the BDS Ansatz 
we would find a similar equation with a slightly different F 

]imd^^^Zn+i{T,e)Rn+i,k{r,e) = dr' F{t,t' ,e){tP{T')Wn,k). (6.3) 

The dependence on e of the OPE coefficient F{T,T',e) is governed, in the e — )• limit, by a 
renormalization group equation which we describe in subsection 6.3. The essential conclusion 
is that the total integral dr does not renormalize, e.g., is e-independent in the limit. 
Thus, for the total integral, which enters Eq. (6.1), F — )• F{a) depends only on the coupling, 
allowing us to write 

1 - f rT=oo 

-^Q{Wn,k) = J Jim J d^^'^Zn+iir, e)i?„+i,fc(r, e) + cyclic. (6.4) 

An important subtlety at this point is that the integral Jq°° dr is singular due to end- 
point divergences. As discussed in the next subsection, the integral is always at most single- 
logarithmic divergent, reflecting the behavior expected for Q of the logarithm of an amplitude. 
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What this means is that in a given ultraviolet regularization scheme the result will be well- 
defined, but it may depend on the scheme.^ 

A natural way to remove the scheme dependence is to divide by the BDS Ansatz, e.g. 
push the factor on the left-hand side inside the Q. Since the BDS Ansatz is one-loop 

exact and proportional to a in the exponent, and since Q is first order in derivatives, this 
adds a term 

{Wn,k)Q^^ = -aRn,k J lim J _^ S^Zn+iir, e)R'^lX^ + cyclic. (6.5) 

Adding the two equations Eq. (6.4) and Eq. (6.5) gives Eq. (1.2), up to the yet undermined 
function of the coupling F{a). In the next subsection we will determine that / F{a) = a, 
using the known fact that QRn,k niust be finite and scheme-independent to all loops [25]. 
The point is that some cancelations are required to occur between the two terms. 

There are various points in this derivation which may not be fully rigorous. For instance, 
we have assumed that a supersymmetric regularization of Wilson loop existed, but, as pointed 
out in [54], in the only regulator scheme which has been tried so far it may be necessary to 
add complicated counterterms to define the correct operator at the quantum level, changing 
the explicit form of the Wilson loop. This means that our derivation is not based on any 
explicitly known regulator. On the other hand, the explicit form of the operator was not 
really important for the derivation; we only really used simple physical properties about the 
excitation spectrum of the Wilson loop. Furthermore, in the end, everything is expressed in 
terms of finite and regulator-independent quantities. For these reasons, we believe that this 
derivation is quite robust. 

Following the same steps for Wilson loops in theories with less supersymmetries, one 
would find that the fermion insertion ip no longer appears inside the chiral Wilson loop, so 
the right-hand side of Eq. (6.1) would be a genuine new object. So = 4 SYM is special in 
being the only gauge theory in which all elementary fields circulate in a chiral superconnection. 
Still, one might be able to derive similar equations in other theories by enlarging the class 
of Wilson loops to be considered. From our viewpoint, the hallmark of integrability is not 
the Q equation itself, because as seen in subsection 3.3 it takes one only ever so far, but the 
existence of a similar equation for Q^^\ which is known so far only for planar J\f = 4 SYM. 

From the scattering amplitude perspective, we expect equations similar to Eq. (1.2) to be 
valid at one-loop in other theories as well (paralleling the results of [31]). On the other hand, 
as noted in Introduction, at higher loops it seems quite difficult, at least to the authors, to 
justify the absence of 1 — )• 3, 4, . . . splitting terms in theories that have no local Wilson loop 
dual. 

In the momentum space introduced in subsection 6.3, the anomalous dimension are O(p^) while the 
endpoint divergences produce 1/p in the small p limit. This is why it is still safe to ignore entirely the 
anomalous dimensions in this discussion. 
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6.2 Convergence of the r integral 

Let us consider the second term in Eq. (1.2), the BDS-subtraction term. Near r = it looks 



This is divergent, reflecting the infrared logarithms in the BDS Ansatz. 

Poles in the r-integrand originate from poles in scattering amplitudes, and this divergence 
can be traced to poles l/(n— 2n— lnn+1) and l/(n— lnn+11). Since the amplitude only has 
poles corresponding to physical channels, these are the only two possible poles which could 
contribute. Actually, the second pole blows up in the collinear limit regardless of r (see 
Fig. 2). A constraint from the /c-decreasing collinear limit implies that the coefficient of this 
pole is Rn,k, so this contribution cancels out between this term and a corresponding part in 
Eq. (6.6), even before we take r — >• 0. The first pole, l/(n—2n— lnn+1), blows up when we 
take r — but this correspond to a soft limit of the amplitude, in which its coefficient also 
reduce to Rn,k- The analysis of divergences near r = oo is similar. We conclude that the 
finiteness of the r-integrand, observed empirically in the previous sections, is a general fact 
which will remain true at any value of the coupling thanks to the nice collinear and soft limits 
of BDS-subtracted amplitudes. 

Note that this is only true when the relative coefficient between the two terms is chosen 
as in Eq. (1.2). This is the reason why we believe that /F{a) = a = jTcusp in Eq. (6.4) 
exactly in the coupling. 

6.3 Absence of log e divergences 

In the main text, we are interested in the zero- momentum component (total r integral) of the 
difference between the two terms in Eq. (1.2). In this subsection, we will be be interested in 
the T-dependence of just the fist term. In particular, we wish to understand the loge terms 
which arise in the e — )• limit at fixed r. 

From general field theory one might expect these logarithms to be related to the anoma- 
lous dimensions of local operators insertions on the Wilson loop. In the context of null 
polygonal Wilson loops this was formalized recently [32], and we refer the reader to this ref- 
erence for more background. In the case at hand the key feature, just alluded to, is that the 
only insertions with the correct (bare) twist and quantum numbers are insertions of single 
fermions. The absence of multi-excitation states is a considerable simplification. This means 
that all pertinent operators are labeled by one parameter, the position along the edge, so we 
expect the renormalization group to act as convolution. Actually, the edge has a symmetry 
which is a combination of a longitudinal Lorentz boost and a dilatation leaving the position of 
the two cusps (and the orientation of neighboring segments) unchanged. In our variables this 
is generated by r — )• a^^^r. It follows that the renormalization group equation is diagonalized 



like 




(6.6) 
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in momentum space. The upshot is that we expect 



Um log 



de/e 



log ex (^(p) - 1) + C7(p) 



(6.7) 



where on the left-hand side, the r integral has been performed but not the e contour integral; 
on the right-hand side, the so-called form factor C(p) (which depends on helicity choices) is 
finite as e — 0, and the so-called dispersion relation Eiji) has to match that of an elementary 
fermion excitation of the null edge (equivalent to excitations of the GKP string [55]), known 
exactly to all values of the coupling thanks to integrability [33] (see also Appendix B of [56]). 

The cancelation of log e divergences at zero-momentum is very easy to understand from 
this formula: the energy £^(0) = 1 is protected by Goldstone's theorem, the zero-momentum 
fermion being the Goldstone fermion for the breaking of supersymmetry caused by the Wil- 
son loop background. The condition -£(0) = 1 is also verified within the integrability frame- 
work [33]. This shows that the cancelations observed empirically in sections 4 and 5 are 
general and will hold exactly in the coupling. 

6.4 The fermion dispersion relation 

We now wish to check the prediction for E{p) at finite p. Due to the physical origin of the 
divergences, it should suffice to check this for the collinear limit of six-point amplitudes, the 
dispersion relation being expected to be universal. We let 

(4512) dT 



lim^d2z„+i 1 d^l^Xn+ii^^tr,? :=Qlog 



(4513) r ^ ' 



(6.8) 



4J- cuspji 



and n = 5. 



as in subsection 4.1, where I is the loop order (e.g., the order in a 

Taking the collinear limit of the six-point tree amplitude (1) + (3) + (5) (the rules in 
subsection 4.1 can be useful here) gives 

1 



rtree 



(r) 



r + 1 



We will need the Fourier transform 
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i sinh ^ 



(6.9) 



(6.10) 



where t = e'^ . The specific form of this result is very useful, as it allows us to write the ratio 

dae^"/"°°P(e"), (6.11) 



where C is the rectangle contour 
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This is vahd for —2 < Imp < 0, where the contributions from infinity can be neglected. Now, 
for any i, /^"l°°P(r) is an analytic function of r with branch points at r = 0, —1, and oo. This 
allows the contour to be deformed and expressed in terms of discontinuities on the horizontal 
line Im cj = vr: 



j£-loop(^) 



dae 2 '^Disc 

-oo 

1 dx 



r^-loop / 



-x^Disc 







X 



-loop f 



(-x) + {p -p) 



(6.12) 
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^ ^ . it could be very mterestmg to interpret 

0, U2 = 1 — ^3, U3 < 0) as the imaginary 



where Disc [/"°°p(-x)] 

this discontinuity (in the cross-ratio regime ui 
part of the six-gluon amplitude in a physical Minkowski-signature regime. 

To deal with such integrals, we found useful to think of y as a positive integer and use 
the language of harmonic sums. Recalling the 5-point r integrand Eq. (4.4) and taking its 
discontinuity, we obtain 



jl-loop 



M 



2x 


loge- 


f M 


[{x-1 


+ J 


loge + 


Si 
N 


-sf 


-S2- 



x + 1 



(x-1). 

.2^ 
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(6.13) 



^ , and the + prescription is the usual 



where M[f] := ^x^ f is the Melhn transform, N : 
one, such that ^^^^l--^^^ = for a > 0. For integer the harmonic sums are defined as 

Si = Yln=l ^ ^"^"i Sij_^i^ = 



, and elsewhere by analytic continuation [57] . The 



— TT /6 term follows from a careful treatment of the x near 1 region, but can also be verified 
numerically quite unambiguously using the first form in Eq. (6.11). 

At two-loops, we are interested only in the loge-terms. Conveniently, these can be read 
off from our formula for the differential of the NMHV hexagon, in Appendix B, without doing 
any integration. The point is that the logarithms exclusively arise from terms proportional to 
dlogui. So, all we have to do, is take the expressions in Appendix B, drop all terms except 
those proportional to dlogui, and expand the degree three functions in powers of logni as 
ui — )• 0. Then we use the simple rules loguidlogui — )• 2 log^ e, dlog ui — )• 2 loge. We then 
have to take a discontinuity as a function of r — )• — x. We obtain for the log^ e terms 
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and for the single-logarithmic terms 



I j2-loop / 

-— = loge M 

4 Jtroc ° I 



+ Li2(x)) + l{l + x) log(l - 3;)(31og(l - x) - logx) 

(1-^)+ 



+ ^C(3)) +(A^^-iV) 



loge f-5i,2 - SiS2 -Sf + ^ + l^-^^-'^Si- |c(3) ) + (iV ^ -N). (6.15) 



iV 2 iV 2iV2 2 2 
The OPE prediction concerns the logarithm, so we need to add the combination 

- I = +loge (25,2 - ^3 + 5? - I - ^1 + ^ + y^i + 3C(3)) + (iV ^ -N) 

-log' el^Sl + - I + ^) + (AT ^ -AT). (6.16) 

This is the square of Eq. (6.13), although writing it as harmonic sums was not entirely trivial 
due to cross-terms between the and —N terms. We found that an efficient way to achieve 
this was to match the poles on the negative N axis. The log^ terms cancel in the sum, and 
we obtain for the logarithm 

^ (log/(p))2-'°°P = loge (^Si,2 -S,-SiS2 + l^ + ^C(3)) +{N^ -N) + finite. (6.17) 

This is to be compared with {E{p) — l)loge where for E{p) we use the expansion to 
second order in Tcusp of the "IctrgG fGrmion" dispersion relo-tion in Ec[.(20) from [33]; 

E{p) - 1 = r,,sp (V'+ - V'(i)) - %^ (v-: + 4V-(^- - ^) + 6C(3)) 

= ^rcuspSi + JrLp(Si,2 - 53 - 52^1 + ^ + 2C(3)) + {N^ -N) 

where V+ := ^(^^(1 + f ) + ^'(l - f )) and V- := i(V'(l + f ) - V'(l - f )) m the first line, 
and on the second line we have converted the result to harmonic sums. Perhaps we should 
have said earlier, that the one-loop prediction 25"! was matched by Eq. (6.13). 

The perfect agreement confirms beautifully that our dP'^^Z integral is probing fermion 
excitations on the edges of the Wilson loop, as was expected from the OPE analysis. Second, 
and perhaps more importantly for us, it gives an independent confirmation to two-loop ac- 
curacy (besides the numerical check in subsection 5.1), that the prefactor in Eq. (1.2) has to 

bp r 

"-"C J- cusp- 

7 Two-dimensional kinematics 
7.1 Preliminaries 

It can be useful to consider special kinematic configurations in which scattering ampli- 
tudes/Wilson loops usually simplify. Following [58, 59] at strong coupling, and [60, 61] 
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at weak coupling, we now consider configurations of external momenta/edges of the Wilson 
loop which can be embedded inside a two-dimensional subspace of Minkowski space. This 
reduces the conformal group SU(2,2) to SL(2)xSL(2). Actually, we are interested in super- 
amplitudes, and as will become apparent soon it is very natural to consider a supersymmetric 
reduction to SU(2,2|4) to SL(2|2) xSL(2|2). 

With no loss of generality we take the number of particles/edges to be even, say 2n. Even 
and odd labels are distinguished, 

Z2^-l = (ai+,o,a2+,o,x^,o,x-+,o), = (o,ai-,o,a2-,o,x!~,o,x-"), (7.1) 

where i = 1, . . . , n for even and odd labels. Four-brackets with two odd and two even labels 
factorize, (2i— 1 2j— 1 2A; 2/) := {ij)[kl], and all others vanish. R-invariants which contain a 
generic reference twistor Z^, , two odd and two even labels will appear, and they factorize into 
odd and even parts, [* 2i—l 2j— 1 2k 21] = (* i j)[* k I] where 

_^^^((HM (7 2) 

^ {*^){^J){J*) ^^^ ^ 

and similarly for the even part [* A; /] . These R-invariants satisfy a four-term identity (a 6 c) — 
(abd) + (acd) — (bed) = 0. Finally, cross-ratios are also separated into odd and even sectors 

In this notation, the NMHV tree amplitude (most easily extracted from the CSW form [63]) 



IS 



7.2 Collinear limits 

The BDS-subtracted amplitudes behave simply under single-collinear limits. In two-dimensional 
kinematics, the natural limit instead collapses the length of an edge to zero, the so-called 
triple-soft-collinear limit. The behavior of amplitudes in this limit is well understood at both 
two-loops and at strong coupling [60, 61]; let us recall the main conclusion. First, it is easy 
to see that nothing can diverge in this limit, just by dual conformal symmetry, to any loop 
order: for the hexagon the whole (BDS-subtracted) amplitude is just a constant since all 
cross-ratios are equal to 1. On the other hand, as is especially clear from the Wilson loop 
viewpoint where the limit has an OPE interpretation [32, 56], the limit involves only physics 
localized around one edge and so cannot depend on the number of points. So in general we 
expect a simple multiplicative relation i?2n — ^ f{o)R2n-2 for some function of the coupling. 

Actually, we will need both ^-preserving or fc-decreasing limits, which are not parity 
conjugate to each other in two-dimensional kinematics. This allows for two different constants. 
It will be useful to absorb them by a simple rescaling of the BDS-subtracted amplitudes: 

R2n,k := e("-2)^('')+'=/2('^)ii2„,fc. (7.5) 
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By choosing /i and /2 such that i?6,o = 1 and iie.i = -^e'^Fi then R2n,k R2n-2,k in the 
/^-preserving hmit X~^^ A", and (/ cPxtd'^Xn^2n,k)/{J d'^Xnd'^Xn(.n-'^'^'^)[n-^nl]) 
i?2n-2,fc-i in the /c-decreasing case. It is known that 

/i(a) = -a2^ + 0(a3), f,(a) = -a"^ + a'^-^ + 0{a'). (7.6) 

The strong coupHng hmits are also known: /i(a) = + /2(a) = + 0(1), when 

extract from [62]. The result for fi at weak coupling has been 
known for a while [60], while the two-loop correction to /2 follows from the value of 
the two- loop value of i?6,i in Eq. (5.4), also obtained recently by [35]. It would be nice to 
have a way to calculate these functions at all values of the coupling. Below we will discuss 
the functions i?2n,fc directly. 

7.3 Two-dimensional Yangian equations 

Eq. (1.2) involves a (2n-|-l)-gon which is not very natural in two-dimensional kinematics. 
However, the key ingredient in its derivation was the fact that the collinear limit of the 
(2n-|-l)-gon had a nonzero overlap with the insertion of a fermion on the 2n-gon. This 
suggests that we can get the same information out of the (2re-|-2)-gon. Namely, we write 
down a limit with the right quantum numbers. Take A and a to be even and add particles 
2n+l and 2n+2: 

QaR2n,k = a j d^^^X+^,J d"|lA;+i(i?2„+2,fc+l-i?2n;2,1^2n,fc)+ cyclic 

■= a Xna _lini _ / (K+idK+i) [ d^x;[+i(t?Xn+i)^ (parenthesis) cyclic. (7.7) 



In principle, one might expect a nontrivial function of the coupling multiplying the first term, 
arising as an OPE coefficient similar to g'^/F{a) above, but not in front of the second, it being 
associated with the BDS Ansatz. However that function would have to be independent of 
n, and arguing as follows it is possible to see that using R this function can only be 1. The 
essential point is that some cancelations have to occur between the two terms. 
On R-invariants not depending n+1, the integral gives zero, while in general 

d'^'Xt^,d'^'X~^,{n+lij)[n+lkl]=Qlogf^dlog^^^^ (7.8) 

L -I \^n+l j / 

The limit in Eq. (7.7) does not depend on how X~_^_i approaches A~, but for individual terms 
it may, thus it is useful to choose A~^^ = A~ -|- eX^ supersymmetrically. Then the previous 
equation is valid in all cases, with the substitution [nn] — )• [nl] when k = n or I = n. The 
action on NMHV tree gives 

d'^'xUd'^'X-+,RlT+2,i = X:<ilog(A+_,,A+)Qlog (7.9) 
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This agrees, for the generic term, with Eq.(3.13) of [61] for the BDS amphtude, but we see 
that the measure is logarithmically divergent at the endpoint X'^_^_i = A^. On the other 
hand, thanks to the universal A;-decreasing collinear limits of R, this divergence will cancel 
in the combination Eq. (7.7) as in subsection 6.2, showing that it is the correct combination. 
Identical equations apply when both indices on Q are in the odd (plus) sector, if we take 

The Q^^^ equation is trickier, since the two-dimensional subsector is not closed under the 
action of the naive generator, Q^^\ We hope to address this issue, and make use of both Q 
and Q^^^ equations in restricted kinematics in the future. 

7.4 From tree N^MHV to one-loop NMHV to two-loop MHV 

Let us now check whether it is indeed possible to compute the two-loop amplitude in restricted 
kinematics, starting with tree amplitudes in only restricted kinematics. For the N^MHV tree 
amplitude we start from the CSW representation [63], which has a simple two-dimensional 
limit for the generic term, albeit the bondary terms are a bit complicated. After some 
massaging we managed to obtain the form: 



R2n% = \ E {*ij){*kl){[ij-lj]-[i-l3-lj]){[kl-ll]-[k-ll-ll]) 

i<j<k<l<i 

I E {*ij)i*ji){{[i^-^i]-[j-ii-imji-^i]-[j-^i-^i]) 

- - [i-ll-ll])) . (7.10) 



3 

i<j<l<i 



Computing the Q from edge 2n using Eq. (7.7) and Eq. (7.10) is a bit tedious, possibly 
because this form for the trees is not optimally simple. After patient bookkeeping we obtain 
the simple result 

Q^2n:r= E ^og^Qlogj^(jkl)i[jl-ll]-[j-ll-ll]) + cyclic. (7.11) 

l<j<k<l<n ' J 

To integrate the Q we need to complete its arguments into cross-ratios. A successful strategy 
is to decompose the logarithms as a difference of two terms and collect the terms. Then for 
the generic term we immediately get cross-ratios, but we also get some boundary terms 

QRlnT = E log{ik){QlogU-_,,_,){jkl){[jl-U] - [i-l/-l/]) - E ^Og{ik){ijk) 

i<j<k<l<i i<j<k<i 

Qlogip^([ii-li] - [i-l j -I j]) + Q log ^-^Tj— ^{[kj -I j] - [k-lj-lj]) 

To complete the arguments of Q in the boundary terms one can in principle follow the 
systematic strategy of subsection 3.4. However in this case this is not necessary, as close 
inspection reveals that adding zero in the form Q log pjj [i j k] = Q log ^ {[k j—1 — j—1 j] + 
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[ij — lk]), antisymmetrized in and {k,k—l), will do the trick. Thus, without much 

attempt at simplifying the result, we obtain the formula 



^2n:r = E umji-^i] - [j-ii-ii])iog{ik)iogu7_^,_, - J2 ^ogmiijk) 

i<j<k<l<i i<j<k<i 

+ ((logt.7.^._, Jii-lfc] - o i-l)) -{k^ k-l))] . (7.12) 



This is the general one-loop NMHV amplitude in two-dimensional kinematics.^ We remark 
that the formula naively contains ill-defined log[jj] terms from k = j+1 in the last sum, 
however these cancel within the square bracket making the formula well-defined. Furthermore 
the formula respects conformal invariance, although it is not expressed explicitly in terms of 
cross-rations. We have verified that this formula agrees with the two-dimension limit of the 
box expansion for 2n = 8, 10, 12. 

To get the two-loop MHV amplitudes from this using the Q equation is much simpler, 
as Eq. (7.8) turns out to reduce to the substitution {*n+lj)[*n+lk] — )• log |^ log[nA;] in this 
case; this yields 

^2nfi^ = - Yl ^°s{ik)log{jl)logur_^j^_-^loguJ_-^^_^-2 ^ log(ij) log(jA:) 

i<j<k<l<i i<j<k<i 

x{loguJ^^^^_^ log^fc-i,i,i-i,i + log^ili,i-i ^oguJ_^ f^ f^_-^ . + lognr_^_^_^ logu^;^.^_i fe_i) 
- 5^ log2(ij) logur_^^^_, log(l - nr ^ J. (7.13) 

i<j<i 

where in the parenthesis the second and third terms are permutations of the first one. We 
have verified that this result is equivalent to the forms given in [56, 61], albeit its form is not 
immediately as elegant. 

Something unusual has to be said. We found that obtaining a reasonably intelligible form 
of the N^MHV tree amplitude in restricted kinematics was by far the most time-consuming 
part of this computation. Obtaining the NMHV one-loop result from it, involved some patient 
bookkeeping, but no special difficulty. Finally, obtaining the two-loop MHV formula was by 
far the simplest step. While we hope that simpler expressions for the tree and one-loop 
amplitudes will be obtained in the future, our main objective here was to confirm that it was 
possible to obtain these results using only input from two-dimensional amplitudes. 



8 Conclusion 

We have proposed an all-loop equation for the Q symmetry acting on BDS-subtracted planar 
S-matrix of = 4 SYM, as well as its parity-conjugate for Q^^^ symmetry, which, interpreted 

'^The first arXiv submission of tfiis paper contained a typograpfiic mistake, witfi [k o fc— 1) appearing as 
(j -o j— 1) in tlie last line. 
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as quantum-corrected symmetry generators, amounts to exact Yangian symmetry. In princi- 
ple, these equations can be used to determine all-loop S-matrix uniquely, at at any value of 
the coupling. As a perturbative study, we have applied the Q equation to reproduce results 
for MHV n-gon, NMHV hexagon at two loops, and fix the undetermined coefficients in an 
Ansatz for three-loop MHV hexagon. The equations relate dual Wilson loops with fermion 
insertions along edges to higher-point ones in the collinear limit, and we have outlined a 
derivation in the spirit of OPE. In particular, we have reproduced the fermion dispersion 
relations to the second order of Tcuspi as a strong consistency check for the proposal. 

Our main finding is that scattering amplitudes are precisely Yangian invariant, provided 
that the correct quantum mechanical expressions for the generator are used, and are fully 
determined by that. 

It should be possible to extend our results in the more loops/more legs directions, with a 
manageable amount of time. For instance, higher-point three-loop MHV amplitudes could be 
analyzed. Using parity the N^MHV heptagon is equivalent to the NMHV heptagon (whose 
symbol we have been able to compute at two loops, though not its function), and from it 
the three-loop NMHV hexagon could be obtained. The two-loop N^MHV octagon should 
be uniquely determined using parity symmetry together with the Q equation; from this one 
could obtain all three-loop heptagons and four-loop hexagons. With, perhaps, the help of 
a big enough computer, we foresee no essential difficulty in obtaining the symbol of these 
objects. Without doing detailed computations, it might also be possible to make general 
qualitative statements about these objects. 

In addition to two- and three-loop results, we have also obtained certain all-loop predic- 
tions for prefactors multiplying degree-(2£— 1) functions for the Q of £-loop NMHV hexagon 
and heptagon, which give a list of the last entry of their symbols. Based on classifications of 
residues in the Grassmannian G(2,n), it should be straightforward to generalize such predic- 
tions to NMHV n-gon, in analog of that of MHV n-gon. More importantly, by considering 
the Q^^^ equation, we hope to carry such analysis for N'^MHV cases, which would provide 
useful information on the all-loop structure of general amplitudes. 

Corresponding functions, as opposed to symbols, could also be obtained provided one is 
able to carry out the one-dimensional integrals which appear, but this might require some 
new ideas. Clearly it would be important to better understand results which are already 
known, for instance the differential of the general MHV n-gon [24] and the three-loop MHV 
hexagon. Indeed the undetermined coefficients in the Ansatz of [36] at three-loops are now 
completely fixed, and it would be fascinating to see the corresponding function. 

It would also be fascinating to study our equations at strong coupling, paralleling the 
the strong coupling application of OPE. It is consistent with all which is presently known at 
strong coupling, in that it express QlogRnfl as \/A times a ratio function, which is known 
to be of order 1. However, to test it in a nontrivial way, one would have to obtain the ratio 
function, which corresponds to a one-loop computation in the string theory. 

Amplitudes restricted to a two-dimensional subspace of four-dimensions seem especially 
promising. There is the intriguing possibility, that the main complexities of loop amplitudes 
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in these kinematics might aheady be present in some way inside the tree amphtudes, where 
the combinatorics have been worked out in [64]. With some effort, one should be able to 
obtain all three-loop MHV amplitudes from some nice form of tree N^MHV. The simplest 
case would be to derive and fix the coefficients of the recent conjecture for the octagon [65] 
from tree N^MHV tetradecagon. 
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A Taking differentials of one-dimensional integrals 

Consider an integral of the form 

/ dlog{x + a)F{x,Ui), (A.l) 

or, more generally, a linear combination of such integrals such that the total integrand con- 
verges both at zero and infinity. 

Then its differential is the sum of the following terms: 

• A term 

- F{x=0,Ui)dloga (A. 2) 

• For each term of the form Gj{x, Ui)dlog{x + xj) in the differential of F, a term 

+ {dlog{a - Xj)) / (dlog— )Gj{x,Ui) (A.3) 

JO ^ + 

• For each term of the form Hj (x, Uj)dlog f in the differential of F, where / is independent 
of X, a term 

POO 

+ {dlogf) / d\og{x-a)Go{x,Ui). (A.4) 
Jo 

The proof is left to the reader; it is more or less integration by parts. Sometimes it may 
happen that some Xj = 0, in which case some intermediate expressions will be ill-defined. 
This can be dealt with efficiently by moving the boundary to x = e, the e-dependence then 
canceling at the end provided the integral is convergent. 

This algorithm for computing derivatives can be used to efficiently (and automatically) 
compute the symbol of any one-dimensional integral; the result depends only on the symbol 
of F, not on its functional representative. 

B Special functions for MHV and NMHV hexagons 

The six-dimensional massless hexagon integral is [45, 46]: 

II"" = --J^ - \J + 2Y,L-{x+Ui,x-Ui), 

i=l 

3 

J = Y.{h{x+Ui)-h{x-Ui)), 

i=l 

+ _ (1245) (3461) (2356) _ (1245)(3461)(2356) 

^ ~ (1234) (3456) (5612)' ^ ~ (6123)(4561)(2345) ' ^ ' 
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This expression is valid for all Euclidean kinematics, ui,U2, U3 real and positive. There appear 
the special combinations introduced in [34] 

LUx^.X-) = ^) ' + ^ i^log(x+X-)- {en-m{x^)+in-m{x~)) , 

L„(X+,X-) = ^ i_A_log(x+X-r {in-m{x-^)-in-mix~)) , (B.2) 

where 

in{x) = -{Unix) - (-l)"Li„(l/x)). (B.3) 

Regarding these special functions, it is important to note that in our applications the ar- 
guments and x~ are the two roots of a quadratic equation with real coefficients. When 
the two roots are complex, the standard branch of the polylogarithms is to be used (the one 
defined on the complex plane minus the line [1, 00)). When the two roots become real, one is 
instructed to add opposite infinitesimal quantities x"*" — )• x"*" + ie, x~ — )• x~ — ie; the hexagon 
function will not depend on the sign of e. 

The colhnear hmit of the N^MHV one-loop amplitude gives rise to following seven inte- 
grals related to the two- loop NMHV 6-point amplitudes: 



QijNMHV ^ (6)Qlog||/i + ((1) - (2) + (4) - (5))Qlog||/2 + ((2) - (4))Qlog||/3 
+ ((6)Q log II + ((5) - (4))0 log h + ((2) + (4))Q log ||/5 

+(5)Q log + (3)0 log M/,, (B.4) 

(63) (63) 

where (1) notates the R-invariant [23456] and 

1 vr^ 
/i = 2^6^ + Li3(l - U2) + Li3(l - I/U2) - Li3(l - ni) - Li3(l - l/m) + — logU3 

+ -(Ll2(l - Ml) +Ll2(l - U2) + Ll2(l - Us)) log \ 

2 U2U3 D 

1 \ 3 1 

Li2(l - U2) + -Li2(l - -^3) j logn3 - -logui log^ "3 - ^ log?i2 log^ U3, 

f2 = Li3(l -U3) (^12(1 -U3) + ^ log'U3logUin2 - logUllogU2) logU3, 

1 o TT^ 

h = -Li3(l - I/U3) + -log'^us + — logn3, 

D 

Ja = ^Li2(l - ui) + Li2(l - U2) + Li2(l - ^3) + logui logn3 - y) logn3, 
/5 = fLi2(l -n3) - log— , /e = log^U3log— , A = logU2logn3log — . (B.5) 
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Expressing the six-gluon NMHV amplitude as 

i?6,i = \ ([(1) + (4)]F3 + [(2) + (5)1^1 + [(3) + {Q)]V2 

+ [(1) - (4)]F3 + [(5) - (2)]Fi + [(3) - (6)]V2) , (B.6) 
we obtain from this, as explained in the main text, the differential 

dVz = -\ird\og ^ + {dV^)id\og . + (dV3)2dlog + (dV3)3C^log 

U2U3 

where 

_ (1235)(2346)(1456) _ (2345)(1356)(1246) _ (1345) (2456) (1236) 

~ (1234) (2456) (1356)' ~ (1345) (2346) (1256) ' ~ (1235) (3456) (1246) ' 

1 Q 

(ciFa)i = 2Lia(l - U2) + Lia(l - 1/1x2) + 2Lia(l - «a) + Lia(l - 1/ua) - - log'^(u2Ua) 
+Li2(l - ui) \0gu2U3 + Li2(l - U2) log — + Li2(l - U3) log — 

Ul Ul 

+- log Ul (log^ U2'ua + 2 log U2 log ua) + — log , 

{dV^)2 = Li3(l - ni) + Li3(l - l/m) + Li3(l - U2) + 2Li3(l - I/U2) - 3Li3(l - Ug) - 2Li3(l - IM) 

+ ^Li2(l -■Ul)log-^ + ^Li2(l -U2)l0g-^ + ^Li2(l -113) log ^^^^ + ^logU2'U3 
2 UlU^ 2 ^21*3 2 o 



1,3 1,3 1 , ,2 1 , , 1 '"2 '''"^ , U1U2 
+ o log" "l^-^S - o log ^2-7; log U1U2 log ^1^3 + - log 1/1 log «3 log ^ - — log , 

3 3 2 2 ^3 2 ua 

(c?^3)3 = ((^1^)2 with ll2 -H- ^3, 

((^^3)4 = - ^Li2(l - Ul) + Li2(l - U2) + Li2(l - 1*3) + log ui log n2U3 - log ^2 log 1*3 - logui 
and 

d% = llfdlog '''^^"'''^ + (dV3)idlogyi + (dV3)2dlogy2y3 + (^^^3)3^^ log -, 
2 (1 - ^2)1x3 y3 

(dF3)i = Li3(l - 1/1x2) - Li3(l - 1/1x3) + f Li2(l - 1x1) - ^"j log ^ - J log^ ^ 

\ d / 1X3 D II3 

(c?V3)2 = -Li3(l - U2) + Li3(l - 1x3) + ^Li2(l - U2) log - ^Li2(l - 1x3) log 

2 U3 2 1X2 

1 / . ,^ , 1 , , 7r^\ 1(2 



+ - |^Li2(l - IXl) + - log IXl log 1X21X3 - yj log 

{dVri)z = Us{l - 1x1) +Li3(l - 1/1x1) - f ^Li2(l - ui) + ^log^ — + ^log^iti + ^ ) log 1x1. 



2 ' ^ 4 ° 1X3 6 ° 6 

These differentials are integrable, e.g., d"^ = 0. The other ones are obtained by cyclic symme- 
try. 
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